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Calculations using staggered quarks augmented with a root of the fermion determinant to reduce 
doubling give a qualitatively incorrect behavior in the small quark mass region. Attempts to 
circumvent this problem for the continuum limit involve an unproven combination of unphysical 
states, a loss of unitarity, and a rather peculiar non-commutation of limits. 
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Figure 1: The phase diagram for one flavor quark-gluon dynamics as a function of m\ and defined in the 
text. A first order phase transition line ending at a second order phase transition appears along the negative 
mass axis. The transition represents a spontaneous breaking of CP symmetry. A finite gap separates the 
transition from vanishing quark mass. 

A popular approximation [0, ||] in lattice gauge theory arises from the simplicity of the stag- 
gered fermion formulation @ ||, §]. With only one Dirac component on each site, the large ma- 
trix inversions involved with conventional algorithms are substantially faster than with other ap- 
proaches. However the approach and its generalizations are based on a discretization method that 
inherently requires a multiple of four fundamental fermions. The reasons for this are related to the 
cancellation of chiral anomalies. To apply the technique to the physical situation of two light and 
one intermediate mass quark requires an extrapolation down in the number of fermions. 

As usually implemented, the approach involves taking a root of the fermion determinant inside 
standard hybrid Monte Carlo simulation algorithms. This step has not been justified theoretically. 
Here I argue that at finite lattice spacing this gives a qualitatively incorrect behavior as a func- 
tion of the quark mass. To produce a correct continuum limit requires the removal of unphysical 
singularities which persist at finite volume and in regions of parameter space where there are no 
physical massless particles. Even in regions where the fermion determinant is positive definite, the 
procedure imposes non-trivial Ward identities incorporating unphysical states. 

The basic problem arises from the exact chiral symmetry of staggered quarks. This symmetry 
appears for each flavor independently and survives the rooting procedure. Due to well known 
anomalies, this symmetry can exceed that allowed in the target theory. 

The issues are sharpest in one flavor QCD. This theory has a rather interesting phase diagram 
when considered as a function of the quark mass [@, fTl ||]. In particular, consider the generalized 
mass term in continuum notation 

mi \jf\l/ + im 5 \j/Y 5 \j/ (1) 

As a function of m\ and m=,, this model has a first order phase transition along a portion of the 
negative m\ axis, as sketched in Fig.(|l|). In many discussions m\ and wis are combined to define a 
"complex" mass m = m\ + im$. 

The nis term introduces an explicit breaking of CP symmetry. When 1115 vanishes, the first order 
transition along the negative mass axis represents a spontaneous breakdown of this symmetry. A 
convenient order parameter is the pseudoscalar expectation 

(V') = i(¥75¥)- (2) 
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This quantity displays a finite jump in sign as one passes through the transition. 

This phase diagram illustrates several important features of one flavor QCD. First, the naive 
classical symmetry under m — > — m is lost, having been broken by the anomaly. Second, this 
transition line does not start at m = but is separated from that point by a finite gap. Indeed, 
physics is expected to be analytic in the mass for a finite region around the origin. The main point 
I am raising is that these features are missing with rooted staggered quarks at finite lattice spacing. 
For the correct physics to be recovered in the continuum would require a complex and implausible 
cancellation of unphysical effects. 

I briefly review the rooted staggered quark approach. Starting with naive lattice fermions, well 
known doubling issues give the fermion spectrum a sixteen fold degeneracy. The staggered ap- 
proach reduces this degeneracy using the observation that for any closed fermion loop, the product 
of gamma matrices associated with that loop is proportional to the identity. If a fermion loop starts 
with only one non-vanishing component on any site, it always involves only a single component on 
each site it passes through. The sixteen doublers immediately factor into four independent pieces. 
This can be formalized by introducing a projection operator on each site 



W^PW (3) 
P = P 2 (4) 
TrP = l (5) 



A particular implementation of this is 

P= I^l + / 7l72 (-l)^+^ + / 7 3 74 (-l)^+^ +75 (_l)^+^+X3+^4^ (6) 

where the x, are the integer coordinates of the lattice site in question. This projection reduces the 
degeneracy from sixteen to four. 

To reduce the doubling further requires something additional, and this is where the rooting 
"trick" comes in. Formally the procedure starts with the path integral and replaces the fermion 
determinant with its fourth root, i.e. \D\ — ► [Dj 1 / 4 . In practice this is achieved by inserting a factor 
of 1 /4 in the fermion force term while traversing a hybrid Monte Carlo trajectory. As D still only 
involves a single fermion component on each site, the process remains quite frugal with computer 
time. Indeed, this is the main motivation for the scheme. 

The problems with this stem from an exact symmetry of the unrooted theory that should not 
be there for the one flavor model. To see this symmetry explicitly in terms of the basic parameters, 
generalize the staggered mass term to 

(mi + iS(x) m 2 ) xjf{x)xif{x) (7) 

where S(x) represents the parity of the site in question, 

S( X ) = (_l)*l+*2+*3+*4_ (8) 

With this extension, the fermion determinant and therefore the path integral are exactly invariant 
under the transformation 

mi^ f cos(d) sin(0) ^ (mi , 

(9) 



mi) \ — sin(0) cos(0)/ \rri2 
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As the mass becomes small, the symmetry requires the appearance of a Goldstone boson or parity 
doubling. For the unrooted theory this is not a problem. The various doublers rotate differently 
under chiral symmetry, and this represents a non-singlet axial symmetry. The Goldstone boson is 
one of the expected pseudoscalar bound states of the four fermion "tastes." 

Is there some connection between the parameter m$ in Eq. ([]]) and the parameter mi in Eq. (^)? 
Naively there would seem to be since the projection operator satisfies 

y 5 P = S(x)P (10) 

However for the unrooted theory the rotation in Eq. @ is a flavored axial symmetry; so, there is 
a hidden flavor mauix in Eq. (^|). For the rooted theory this symmetry is spurious, the theories are 
different, and the question of this paragraph is moot. 

To make the problem with rooting more precise, consider working in finite volume and on 
a lattice with an even number of sites. In that case the determinant of the fermion matrix \D\ is a 
polynomial in the combination m\+m\. This is easily seen from the hopping parameter expansion. 
Furthermore \D\ is real and non-vanishing for real m^+m^ > 0. This comes by looking at the 
eigenvalues of D after rotating away mi in favor of mi. Finally, there is a lower bound on the 
determinant 

\D\>(m\ + m 2 2 ) N ^l 2 (11) 

where N c is the number of colors and V the number of lattice sites. 

These facts show that roots of the determinant are analytic in m\ and mi for at least some finite 
vicinity of any point where m\ + m\> 0. This analyticity allows a continuation from m to — m. For 
example one can follow the path 

mi=mcos(0), m2 = msin(0), < 8 < 71. (12) 

This continuation avoids all singularities in jDl 1 / 4 . 

We conclude that for rooted staggered fermions all correlations at mi are identical to those 
at —mi. Furthermore the Ward identities that follow from -jh-qZ = require a Goldstone boson 
or parity doubling. These properties are in direct contradiction to the known behavior of the one 
flavor theory, which is analytic at m = with the m and — m theories being inequivalent and there 
is no Goldstone boson or parity doubling. While the one flavor theory is somewhat peculiar, these 
problems extend directly to more flavors with non-degenerate quarks. This is because one has an 
independent U(l) symmetry for each flavor, which is one more symmetry than allowed. 

The advocates for staggered quarks [§] have argued that these facts would be irrelevant if the 
diseases were to go away in the limit of vanishing lattice spacing. They argue that the the limits 
a — > and m — > might not commute. They point out that instanton effects can bring in \Jm 2 + a 2 
factors which break the m — ► — m symmetry if a — > first. They do allow that there are extra Ward 
identities, but these could be satisfied by indefinite metric "taste" states. The required Goldstone 
boson would be there, but as one of a larger set of unphysical states. The use of an indefinite metric 
could permit a cancellation of these extra states from all physical processes. 

Is this picture for saving the rooting procedure plausible? The idea of infrared and ultraviolet 
limits not commuting seems a bit strange, but the staggered proponents are asking a lot more. The 
problems persist even at finite volume, where one would not expect any infrared issues with the 
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mass going to zero. Furthermore, for the one flavor theory the physical spectrum does not contain 
any massless particles; so, it is unclear why an infrared regulator should be required anyway 

Note that the chiral behavior at finite cutoff a is explicitly wrong because of the extra Goldstone 
boson; indeed, this invalidates what was once regarded as a virtue of staggered fermions, the fact 
that they have some remnant chiral symmetry. Furthermore, the use of an indefinite metric for 
unphysical states requires a major leap from the canonical lattice ideas of Osterwalder and Seiler 
[JlO|]. Finally one should note that these issues are not present with other fermion regulators, such 
as Wilson [O], domain wall [12], or the overlap [O]. 



In conclusion, at finite lattice spacing rooted staggered quarks give a qualitatively incorrect 
description of one flavor QCD. Attempts to circumvent this require complicated and implausible 
limits, not necessary with other lattice regulators. Without proof, rooted staggered quarks should 
be regarded as an uncontrolled approximation, from which claims of first principles calculations 
are at best premature. 
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